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The process 7 + f— ^n + dis treated by means of three-body integral equa- 
tions employing in their kernel the W-Matrix representation of the subsystem 



> 

(N 

^\j ' amplitudes. As compared to the plane wave (Born) approximation the full 



C^ ' solution of the integral equations, which takes into account the final state 

interaction, shows at low energies a 24% enhancement. The calculations are 
Q ' based on the semirealistic Malfliet-Tjon and the realistic Paris and Bonn B 

C ' potentials. For comparison with earlier calculations we also present results 

:^ • for the Yamaguchi potential. In the low-energy region a remarkable potential 

X 

j^ ■ dependence is observed, which vanishes at higher energies. 
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I. INTRODUCTION 

The photodisintegration of '^H and ^He and the inverse reaction, the radiative capture 
of protons or neutrons by deuterons, have been intensively investigated in the past. Due 
to the fact that the corresponding matrix elements contain both the three-body bound and 
continuum states, these reactions are expected to be a good testing ground for the underlying 
two-body potential. In early calculations by Barbour and Phillips ||I|] it has been shown that 
an exact treatment of the continuum states results in a considerable enhancement of the cross 
section in the peak region. In that work and in the following by Gibson and Lehman |^ only 
simple s-wave interactions of Yamaguchi type have been used, while d-components {j < 1^) 
were incorporated in [^. The role of p-wave contributions in the two-body input with 
respect to the three-body cross section was investigated by Fonseca and Lehman Q, using 
again Yamaguchi terms in the interaction. More recently, calculations based on separable 
representations of the Bonn A and Paris potentials, including their higher partial wave 
contributions, were performed for polarization observables at some specific energies 0-|^. 

In the present work we calculate the differential cross sections at 90° for the Yamaguchi, 
Malfliet-Tjon, Paris, and Bonn B potentials over an energy region from threshold up to 40 
MeV. The potential dependence, the effect of higher partial waves, and the role of meson 
exchange currents, taken into account via Siegert's theorem, are investigated. 

Technically the calculations are based on the Faddeev-type Alt-Grassberger-Sandhas 
(AGS) formalism p adjusted to photonuclear processes JlOl, as done already in 0. The 



separable representation of the subsystem T matrices, relevant in this context , is chosen 
according to the W-Matrix approach ||ll| . This approximation combines high accuracy with 
considerable simplicity |1T2| . 

A more reliable, still fairly simple representation of the two-body input is provided by the 
Ernst-Shakin-Thaler method [|1^. Calculations employing the corresponding PEST, BEST, 
and NEST potentials [p^-[T6[| will be presented in a subsequent publication. 



II. FORMALISM 

The Alt-Grassberger-Sandhas (AGS) equations are well known to go over into effec- 
tive two-body Lippmann-Scliwinger equations |^ when representing the input two-body T- 
operators in separable form. The neutron-deuteron off-shell scattering amplitude T{q,q"), 
thus, is determined by 

T{q, q") = V{q, q") + j d\' V(q, q') Goiq') Tiq' , q"). (1) 

Applying the same technique to the photodisintegration of the triton into neutron and 
deuteron, i.e. to the process '-f + t —>■ n + d, an integral equation of rather similar structure 
is obtained 11,10, 



M{q) = B{q) + jd^^ V(g, q') Goiq') M{q') , (2) 

where M. [q) represents an off-shell extension of the full photodisintegration amplitude 

M{q)= (-)(g;^,|//em|^t). (3) 

In both these equations the kernel is given by the effective neutron-deuteron potential V 
and the corresponding effective free Green function Qq defined in Ref. 0. However, in Eq. 
(I) the inhomogeneity of Eq. (P is replaced by an off-shell extension B{q) of the plane-wave 
(Born) amplitude 

B{q)= (gl(V^dli^eml^t). (4) 

Here |\l'i) and {ipd) denote the triton and deuteron states, \q) the momentum state of the 
neutron relative to the center of mass of the deuteron, and i^em the electromagnetic operator. 
Thus, by this replacement any working program for n-d scattering can immediately be 
applied to the above photoprocess. 

The results presented in this paper are, in fact, obtained by extending recent n-d calcu- 
lations [|r3,0] in this manner. As in these references, the separable part of the W-matrix 



representation [|Il|] of the two-body T matrices is employed. 



T;1,{p,p', E + tO) = J2p'w;jip, k; E)Al{E + ^0)W;^r,ip' , k; E)/' . (5) 

iv 

Here / and /' are the orbital angular momenta, and r] = {s,j;i) stands for the spin, the 
angular momentum j [ with the coupling sequence (/, s)j] and the isospin i of the two-body 
subsystem. With the form factors W^ and the propagators A~ of this representation Eq. 
(0) reads, after partial wave decomposition and antisymmetrization. 



Ml{q,E) = ^Bl{q,E) + J] J dq'q" ^V\f: {c^ ,E)A%{E ~ {q'^)^ M^iq^E) (6) 



b'i',i" 
and its inhomogeneity is given by 



oo 

^B'^{q,E) = V3J2 J dpp^+'{pqbiri\W]ik,p,E - lq'')Go{E + iO)H,^\^ilt) . (7) 



' 

Here, Gq and {pqhlTI\ represent the free three-body Green function and the partial wave 
projection of the three-body plane wave state {q\{p\- The detailed structure of the effective 
potential '"V^^-^ can be found in Ref. ||12[. The label b denotes the set ijiKL) of quantum 



numbers, where K and L are the channel spin of the three nucleons [ with the coupling 
sequence (j, \)K] and the relative angular momentum between the two-body subsystem and 
the third particle, respectively. F is the total angular momentum following from the coupling 
sequence {K, L)r, and / is the total isospin. 

The Born amplitude contains the triton bound state |\l/t), which may be calculated by 
any of the various bound-state methods. Consistently with the present approach we employ 
for this purpose the partial wave projected homogeneous version of Eq. (^, 

oo 

Fii<D =T.Jdq' q" Vf '(g, q\ E^) A}{Et - | q") FP'{q'). (8) 

b'i'l 

Its solutions, the three-body "form-factors" Fl'{q), are related to |\l/() according to by 
1**)= Y.JJdqdpq'p'Go{ET) \{^)qpblVI)W,i{p,k,ET - Iq") /^l^Er - \q') F^iq). (9) 

Here, the summation runs over all two-fragment partitions 7, the variables and quantum 
numbers being understood in the corresponding set of Jacobi coordinates. 



The electromagnetic operator entering Eq. (^ is, at the low energies considered, essen- 
tially a dipole operator. Ignoring meson-exchange currents it is given by 

^em = yYEe.P.>lA(^,¥')- (10) 



With exchange currents it takes, according to Siegert's theorem |jT8[, the form |T9 

Hem = -m/— {Ef - Et) Y.e,nY^^{^,ip) . (11) 

V "J i=i 

Here Ef and Et denote the final and the triton energies, r, the nucleon center-of-mass 

coordinates, Pi the corresponding momenta , Cj the electric charges, and A the polarization 

of the photon. 

The on-shell restricted solutions ^M^ of the integral Eq. (^ yield the photodisintegration 
amplitude via the partial-wave summation 

(") ( qsrus; ^djrrij \ i^cm | *t T'Mr' ) 
= J2 ( J^^^^s I KMk ) ( KMkLMl I TMr ) Ylm, (q) ^M\ {q, E, + f g^) , (12) 

rMrblMKML 

where now, in contrast to Eq. @, the spin and angular momentum quantum numbers of the 
neutron [srUs), deuteron (jnij), triton (r'Afp), and the polarization (A) of the photon are 
explicitly given. Denoting this amplitude by M{q)msmjM^,\ , the cross section is obtained 
in the standard way by 



da 2tt^ q 



Y: E M(g )^^^^m„a . (13) 



Here we have averaged over the initial states and summed over the final states. 

III. RESULTS 

As a first test of our numerical program we performed calculations for the Yamaguchi 
PT| potential in order to compare them with the corresponding results by Gibson and 
Lehman 0. Unfortunately, employing the same sets of parameters as in this reference, we 



read off from the analytically given relations somewhat different values for the two-body 
scattering length, effective range, and deuteron binding energy (Table I). The triton binding 
energies, obtained by solving the homogeneous three-body equations, differ in one case, the 
symmetric Tabakin potential, where we found — 9.72MeV instead of — 9.33MeV (Table III). 
These discrepancies mean that our photodisintegration results cannot be expected to fully 
agree with the ones of Ref . [@] . 

Figure 1 compares our cross sections for triton-photodisintegration (solid line) with the 
calculations by Gibson and Lehman |0| (dashed line). The disagreement, in particular in 
the peak region, is reduced when replacing in the Siegert operator our triton energy by the 
Gibson-Lehman value (short-dashed line). Figure 2 shows the same comparison, but now 
for the ^He photodisintegration. Since there is no disagreement between the ^He binding 
energies, a correspondingly modified curve does not exist. The remaining discrepancies, 
therefore, have to be attributed to numerical uncertainties, which are not unexpected in 
view of the level of accuracy reached in early calculations. Within these limits we consider 
our results for the Yamaguchi potential consistent with the ones of Ref. 0. 

The differential cross sections obtained for the Paris p2|, Bonn B P3|, Malfiiet-Tjon 



(MT I-III) [^, and Yamaguchi potentials are shown in Figure 3. Most remarkable is the 
strong potential dependence in the peak region, which vanishes at higher energies. A further 
observation, which should be relevant in model calculations or when going over to higher 
particle numbers [^, is the proximity of the Paris and Malfliet-Tjon results on the one 
hand, and less closely the Bonn B and Yamaguchi potentials on the other hand. 

In Figure 4 we contrast, for the Paris potential, the solution of the integral equation 
(solid line) with the corresponding plane- wave (Born) approximation (dashed line). It is 
seen that the full solution is enhanced by about 24% at the peak. A similar enhancement 
was observed for simpler interactions already in Ref. pi]. The upper curves are the ones based 



on the Siegert-operator (|TT]), the lower ones correspond to the non-Siegert operator (|T0[). 
There is a factor of two between the Siegert and non-Siegert results, which demonstrates the 
relevance of meson-exchange currents. For a detailed discussion of the same phenomenon in 



case of deuteron and ^He photodisintegration we refer to Refs. [^] and [^, respectively. 
Figures 5 and 6 compare our cross sections for the s- and d-wave projected Paris and 



Bonn B potentials with the data from Refs. p6|-E3]. Figures 7 and 8 show the same with 



inclusion of the subsystem p waves, which leads to a reduction of the peak by 8-10%. Up 
to 25 MeV the best agreement with the data is achieved for the Bonn B potential. In 
view of the experimental errors the relevance of this observation is, of course, somewhat 
questionable. At higher energies the potential dependence vanishes. It is, however, seen 
that the incorporation of the p waves is essential for the remarkable agreement with most 
recent experimental data |2^. Note that for the inverse reaction the relevance of the p wave 
contributions has been pointed out also in ||^-§[. 

Two fragment photodisintegration of the three-nucleon bound states, thus, provides a 
sensitive tool for testing the underlying two-nucleon potentials. A repetition of the corre- 
sponding low-energy measurements with much higher accuracy is strongly suggested by this 
observation. 
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TABLES 
TABLE L Parameters for the Yamaguchi-Potential. The numbers in parenthesis are taken 

from [2]. 















Binding 






Strength 


Inverse range 


Scattering length 


Effective range 


energy 


Interaction 




A (fm-3) 


P (fm-i) 


a (fm) 


ro (fm) 


(MeV) 


n-p triplet 


I 


0.3815 


1.406 


5.433 

(5.423) 


1.761 

(1.761) 


-2.203 

(-2.225) 




II 


0.220 


1.15 


5.806 

(5.68) 


2.088 
(2.09) 


-2.082 
(-2.225) 


n-p singlet 


I 


0.1445 


1.153 


-23.196 

(-23.715) 


2.732 

(2.74) 






II 


0.148 


1.15 


-42.217 
(-21.25) 


2.680 

(2.74) 




p-p singlet 




0.1534 


1.223 


-7.853 
(-7.823) 


2.794 
(2.794) 




n-n singlet 




0.1323 


1.130 


-16.851 

(-17.0) 


2.841 

(2.84) 




n-p triplet 


III 


-0.220 


1.14525 


5.666 


2.081 


2.226 


n-n singlet 


III 


-0.148 


1.16225 


-22.998 


2.710 





TABLE II. Parameters for the Malfliet-Tjon-Potential. 



a 



A^[fm- 



A?[fi 



m 



r^ [fm] 



r^[fm] 



-19.5719 

-23.8775 



109.605 



0.643087 



0.321543 



TABLE III. Triton Binding energies obtained with the Yamaguchi- Potentials. The numbers in 
parenthesis are taken from [2]. 



Wave function 



Interaction set 



Binding energy 
(MeV) 



Symmetric Tabakin 



Tabakin 
Charge dependent 

Charge dependent 

(adjusted) 
Charge dependent 

3He 

(adjusted) 



Average of n — p 
triplet and singlet 

sets II 

n — p triplet and 

singlet sets II 

n — p triplet I 

n — p singlet I 

n — n singlet 

n — p triplet 

A = 0.3608 fm-3 

n — p singlet I 

n — n singlet 

n — p triplet I 

A = 0.3589 fm-3 

n — p singlet I 

p — p singlet 



9.72 
(9.33) 

10.11 

(10.1) 

10.34 

(10.34) 

8.49 
(8.49) 



7.72 
(7.72) 



3H 

(present) 



n — p singlet III 
n — n singlet III 



9.968 
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FIGURES 

FIG. 1. Triton photodisintegration cross section for the Yamaguchi potential compared with 

[2]. 

FIG. 2. Same as Fig. 1 but for ^He. 

FIG. 3. Triton photodisintegration cross sections for the Yamaguchi, Malfliet-Tjon, Paris 
and Bonn B potentials. 

FIG. 4. Triton photodisintegration cross sections for the Paris potential {j < 1^) with and 
without Siegert's theorem. 

FIG. 5. Cross sections for the Paris potential (j < 1^) compared with the data from Refs. 

[23-25]. 

FIG. 6. Same as Fig. 5 for the Bonn B potential. 

FIG. 7. Cross sections for the Paris potential with incorporation of the p waves {j < 1) 
compared with the data from Refs. [23-25]. 

FIG. 8. Same as Fig. 7 for the Bonn B potential. 
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